rt<-read.table("tmc.TXT",head=FALSE)
radi<-rt[[2]]
texc<-rt[[3]]
nvol<-rt[[4]]
ncol<-rt[[5]]
sthm<-rt[[6]]
snth<-rt[[7]]
sthd<-rt[[8]]
mlte<-rt[[9]]
mvir<-rt[[10]]
mjns<-rt[[11]]   

lgradi<-log(rt[[2]])
lgtexc<-log(rt[[3]])
lgnvol<-log(rt[[4]])
lgncol<-log(rt[[5]])
lgsthm<-log(rt[[6]])
lgsnth<-log(rt[[7]])
lgsthd<-log(rt[[8]])
lgmlte<-log(rt[[9]])
lgmvir<-log(rt[[10]])
lgmjns<-log(rt[[11]])  

ty=lgmvir
tx=lgmlte
cor.test(ty, tx, method="pearson")